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Abstract 

Background: Williams syndrome transcription factor (WSTF) is a multifaceted protein that is involved in several 
nuclear processes, including replication, transcription, and the DNA damage response. WSTF participates in a 
chromatin-remodeling complex with the ISWI ATPase, SNF2H, and is thought to contribute to the maintenance of 
heterochromatin, including at the human inactive X chromosome (Xi). WSTF is encoded by BAZIB, and is one of 
twenty-eight genes that are hemizygously deleted in the genetic disorder Williams-Beuren syndrome (WBS). 

Results: To explore the function of WSTF, we performed zinc finger nuclease-assisted targeting of the BAZIB gene 
and isolated several independent knockout clones in human cells. Our results show that, while heterochromatin at 
the Xi is unaltered, new inappropriate areas of heterochromatin spontaneously form and resolve throughout the 
nucleus, appearing as large DAPI-dense staining blocks, defined by histone H3 lysine-9 trimethylation and associ- 
ation of the proteins heterochromatin protein 1 and structural maintenance of chromosomes flexible hinge domain 
containing 1. In three independent mutants, the expression of a large number of genes were impacted, both up 
and down, by WSTF loss. 

Conclusions: Given the inappropriate appearance of regions of heterochromatin in BAZIB knockout cells, it is 
evident that WSTF performs a critical role in maintaining chromatin and transcriptional states, a property that is 
likely compromised by WSTF haploinsufficiency in WBS patients. 



Background 

X chromosome inactivation (XCI) is the process whereby 
females balance the levels of X-linked gene expression 
with males [1], and is an archetypal example of epigenetic 
regulation. Much is known about the early stages of XCI, 
including counting, choice, initiation, and spreading of the 
inactivation signal [2], but our understanding of how chro- 
matin of the chosen inactive X chromosome (Xi) and gene 
silencing is faithfully maintained throughout subsequent 
somatic cell divisions is less well understood. 

The human Xi is primarily composed of two spatially 
distinct types of heterochromatin [3]; one is character- 
ized by histone H3 trimethylated at lysine 9 (H3K9me3) 
[4,5] and association of heterochromatin protein 1 
(HPl) [6], whereas the other is defined by histone H3 
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trimethylated at lysine 27 (H3K27me3) [7,8], elevated 
levels of the histone variant macroH2A [9-11], and associ- 
ation of the X-inactive specific transcript (XIST) [12-14], 
We have previously shown that the WSTF-ISWI chroma- 
tin remodeling complex (WICH) transiently associates 
with the human Xi as the chromosome is undergoing 
DNA replication [15] and, therefore, is a candidate for 
maintaining this chromatin organization. Current models 
suggest that Williams syndrome transcription factor 
(WSTF), one of two subunits in the WICH chromatin re- 
modeling complex [16], assists in reforming chromatin 
states of the parental cell post DNA replication [17], Con- 
sistent with this model, depletion of WSTF protein levels 
by RNA-interference (RNAi) resulted in aberrant hetero- 
chromatin formation [18], 

WSTF is encoded by the gene BAZIB, which is lo- 
cated at 7qll,23 [GenBank: NG027679] [19], It is a 
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multifaceted protein with ties to several cellular pro- 
cesses, including transcription [20], replication [16], and 
the DNA damage response [21]. In addition to the 
WICH complex, WSTF has been identified as a partici- 
pant in the B-WICH complex, which is associated with 
transcription of the ribosomal genes [22]. WSTF is also 
linked to the DNA damage response pathway, and was 
recently shown to possess tyrosine kinase activity, phos- 
phorylating Tyrl42 of H2A.X [21]. The significance of 
this covalent histone modification sits at the juncture 
between deciding on DNA repair or apoptosis [23] based 
on whether the modification is removed to allow access 
to the MDCl complex or retained resulting in activation 
of the JNK pathway [24]. 

Interestingly, WSTF is deficient in Williams-Beuren syn- 
drome (WBS; OMIM#194050), an autosomal dominant 
genetic disorder that affects approximately 1 in 7,500 indi- 
viduals [25]. The disorder is characterized by a 1.5-1.8 
megabase (Mb) deletion from one copy of chromosome 
7qll.23 [26]. The deletion arises de novo, due to low-copy 
repeat elements flanking either side of the break points 
[27], and comprises approximately twenty-eight genes 
[28-31]. WBS patients experience an array of symptoms 
including cardiac defects, cognitive impairment, hypercal- 
caemia, growth deficiencies, and a distinct craniofacial 
phenotype [32-34]. Stemming from studies performed in 
mouse and cell lines, WSTF has been implicated as a key 
contributor to many of the phenotypes seen in WBS pa- 
tients, despite the number of other genes that are also de- 
ficient [16,35,36]. It is unclear, however, which specific 
gene deficiencies contribute to each symptom; a problem 
common to all contiguous gene deletion disorders. 

The aim of this study was to investigate the impact of 
WSTF loss on the maintenance of chromatin at the hu- 
man Xi, and more broadly, to develop a better under- 
standing of how the loss of this multifaceted protein 
affects general cellular processes. Towards this goal, we 
disrupted the BAZIB gene using zinc-finger nuclease- 
assisted gene targeting in the cell line hTERT-RPEl 
(RPEl)— a telomerase immortalized [37], diploid female 
cell line— and isolated several independent clones that 
lack the WSTF protein. Findings from this study under- 
score the necessity of determining the role of WSTF, 
particularly in regards to chromatin maintenance. 

Results 

Disruption of the BAZIB gene using zinc finger nuclease 
technology 

The development of zinc finger nuclease technology 
provides opportunities for genomic editing at precise lo- 
cations in human cells [38]. Here, we introduced into 
RPEl cells a pair of zinc finger nucleases (ZFNs) de- 
signed to exon-7 of the BAZIB gene, along with an engi- 
neered repair template (Figure la). The repair template 



consisted of a central promoterless neomycin open read- 
ing frame (ORE) that was preceded by a splice acceptor 
and internal ribosome entry site, and followed by a poly- 
adenylation signal [39]. This central cassette was flanked 
by homology arms (HA) of approximately 600-800 bp 
that were derived from sequences in the vicinity of 
exon-7 (Figure la). The general strategy is that upon 
generation of a double-stranded break by the ZFNs, the 
repair template is used though homology mediated re- 
pair to patch the damage, introducing the neomycin 
ORE in place of exon-7 in the process. This synthetic 
exon promoter trap (SEPT) is incorporated into the tar- 
geted BAZIB transcript, prematurely truncating the 
message through poly-adenylation and providing neomy- 
cin resistance, driven by the genes own expression. 
Transfected cells were placed under neomycin selection, 
and after two weeks clones emerged. 

Screening for correct targeting was achieved by per- 
forming PGR for both the left and right sides of the tar- 
geted region with a primer within the SEPT cassette and 
a primer located outside of the HAs at the genomic 
locus (Figure lb). We screened a total of 119 single cell 
clones by PGR. We identified 13 clones that were cor- 
rectly targeted at one or both alleles of the BAZIB gene. 
To screen for targeting at both alleles, a second PGR 
was performed with primers flanking the ZEN target site 
(Figure la). If a clone was targeted at both alleles, this 
sequence would no longer be present and no product 
would be detected. Three clones were determined to be 
successfully targeted at both alleles. Double-strand 
breaks generated by ZFNs can be repaired by either 
homology mediated repair or non-homologous end join- 
ing (NHEJ). NHEJ is not a precise repair pathway, fre- 
quently resulting in introduction of small insertions or 
deletions at the target site [40], which can disrupt the 
reading frame of the target gene. For those clones that 
were targeted at only one allele, we screened for NHEJ 
at the second allele by sequencing the PGR product. We 
identified three clones that were targeted at one allele 
and have the second allele disrupted by NHEJ (Figure Ic), 
and several clones that were targeted at one allele but 
wild-type at the second allele. This latter group of het- 
erozygous mutant clones are useful models to investigate 
the contribution of BAZIB hemizygosity in WBS in the 
absence of the other >20 genes that are typically lost 
alongside BAZIB in this disease, all of which potentially 
contribute to the phenotype. The clones used in this 
study and the nature of their mutation/s is indicated in 
Table 1. 

Quantitative real time PGR (qRT-PGR) using primers 
that spanned exons-2 and -3 indicated that a clone tar- 
geted at only one allele (Glone A6), showed BAZIB 
mRNA at 47.6% of parental cells. Glone Ml, targeted at 
one allele and mutated at the second by NHEJ, expressed 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 Generation of BAZ1B Icnoclcout cells, (a) Scliematic representation of tine BAZIB gene. Exons: numbered vertical bars. Introns: 
connecting lines. Exon 7 is bracketed by the blue lines, and expanded beneath. The targeting construct is shown beneath, indicating the location 
of the neomycin resistance gene and the extent of the homology arms (regions between the vertical dotted lines). Left and right side screening 
PCRs are indicated below by inward-facing arrowheads. The NHEJ PGR is shown above Exon 7, and the ZEN cutting site is indicated by the white 
arrowheads within the exon. (b) Example of an ethidium bromide stained gel showing a positive clone by PGR. (c) DNA sequence alignment of 
parental RPEl with NHEJ disrupted clones. The cutting site for the ZFNs is highlighted in red. Deletions in NHEJ clones are indicated by the red 
gaps; insertions by the blue triangle, (d) qRT-PGR analysis of BAZIB showing levels in a heterozygous mutant {BAZIB +/-), a clone with one allele 
targeted and one disrupted by NHEJ {BAZIB -/NHEJ) and a clone targeted at both alleles {BAZIB -/-) relative to parental RPEl. (e) Indirect im- 
munofluorescence showing the distribution of WSTF in RPEl, a BAZIB +/-, BAZIB -/NHEJ and a BAZIB -I- clone. WSTF is shown in green for the 
Bethyl Labs anti-WSTF antibody and red for the Gell Signaling anti-WSTF antibody. The white arrowheads indicate the location of the Xi. The nu- 
cleus is counterstained with DAPI (blue), (f) Western blot for WSTF with the Gell Signaling antibody in RPEl, a BAZIB +/-, and various independ- 
ent BAZ]B knockout mutants. Actin is shown immediately below, (g) Schematic map of WSTF showing the relative location of functional 
domains, and the truncated proteins possible in the different mutant types, (h) Western blot using the Bethyl Labs anti-WSTF antibody as 
described for part-f above. 



BAZIB mRNA at 183%, whereas clone D5, that is tar- 
geted at both alleles, showed 33.7% (Figure Id). Expression 
of BAZIB is expected as exons 1-6 are uninterrupted and 
isolation of clones required BAZIB expression to trap the 
SEPT cassette and express neomycin resistance. The re- 
duced levels of mRNA likely represent destabilization of 
the BAZIB mRNA. 

We then employed indirect immunofluorescence (IF) 
to investigate the distribution of WSTF in the BAZIB 
mutants. Here, we used two independent antibodies 
raised to different regions of the WSTF protein; one to 
the first 50 amino acids, and the other to a peptide at 
the carboxy terminus. Results indicate that WSTF distri- 
bution in the BAZIB heterozygous mutant was very 
similar to that of the parental cells, whereas WSTF 
staining was absent in all double-targeted clones sug- 
gesting a complete knockout of the gene (Figure le). To 
confirm loss of WSTF, Western blotting was performed 
on protein extracts from the mutants. Consistent with 
the IF analysis, protein levels were reduced in the het- 
erozygous mutant and absent from all independent 
homozygous mutant clones (Figure If). By targeting 
exon-7, the ORE contained within the first 6 exons could 
potentially generate a truncated protein of 87.3 kDa 
(from NHEJ targeted alleles, amino acids 1-761) or 
34.5 kDa (from double targeted alleles, amino acids 1- 

Table 1 BAZIB mutations in clones used in this study 



Clone 




BAZ1B Allele 1 


BAZ1 B Allele 2 


BAZ1B-A6 {+/-) 


Promoter trap targeted 


Wild type 


BAZ1B-D5 { 


-/-) 


Promoter trap targeted 


Promoter trap targeted 


BAZlB-f3 { 


-/-) 


Promoter trap targeted 


Promoter trap targeted 


BAZlB-m { 


-/-) 


Promoter trap targeted 


16 bp deletion 


BAZ1B-N6 { 


-/-) 


Promoter trap targeted 


2 bp deletion 


BAZlB-04 { 


-/-) 


Promoter trap targeted 


4 bp insertion 



Table indicates the nature of the mutations generated at the two 
chromosome 7 BAZIB alleles for the clones used in this study. Next to each 
clone name, (+/-) indicates heterozygous mutant whereas (-/-) indicates 
homozygous mutant. 



297) (Figure Ig). The NHEJ form would contain the 
WAC and DDT domains only, whereas the double tar- 
geted form of the protein would only contain the WAC 
domain. In order to detect such a protein, we repeated 
the Western analysis using the second antibody that was 
raised to the first 50 amino acids of WSTF. Consistent 
with the first Western, no truncated protein could be 
detected (Figure Ih). 

Loss of WSTF does not impact chromatin of the human Xi 

At interphase, the Xi can be readily detected as a densely 
staining heterochromatin mass, or Barr body, that is typic- 
ally located in a perinucleolar or peripheral location in the 
nucleus [41]. The Xi is painted by XIST RNA [14] and fac- 
ultative heterochromatin markers are readily observed due 
to the chromosome wide silencing, resulting in an un- 
usually large concentration of these features. Previously, 
we have shown that WSTF as part of the WICH complex, 
associates with the Xi territory during S-phase [15] as the 
chromosome is replicated [42]. In order to determine if 
WSTF loss impacts the organization and maintenance of 
Xi chromatin, we assessed a number of Xi markers for 
changes in their distribution. 

First, we performed RNA FISH to investigate whether 
the Xi had retained its association with XIST RNA. 
Here, we found that XIST RNA was enriched in the 
vicinity of the Barr body in BAZIB knockouts in a pat- 
tern indistinguishable from that of parental RPEl cells 
(Figure 2a). Next, we examined other chromatin features 
of the Xi. Typical characteristics of the human Xi in- 
clude an enrichment of H3K9me3 [4,5], H3K27me3 
[7,8], and mH2Al [10], and a noticeable absence of 
H3K4me2 [4,43] and histone acetylation [44], such as 
H3K9Ac. By IF, none of these features appeared per- 
turbed in the BAZIB knockouts, having distribution pat- 
terns indistinguishable from that observed in parental 
RPEl cells (Figure 2b). However, of note, the distribution 
of H3K9me3 throughout the rest of the nucleus ap- 
peared more punctate (Figure 2b, top row). 
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(a) RPE1 BAZ1B -/- 




Figure 2 Investigating the impact of WSTF loss on the human 

Xi. (a) Distribution of XIST RNA in RPEl and BAZIB -I- using a 
direct-labeled RNA FISH probe (green). Nuclei in all images are coun- 
terstained with DAPI (blue), merged images are shown on the right 
of each group of three images, and white arrowheads indicate the 
location of the Xi. (b) Examples of the distribution of markers that 
are elevated at the Xi territory including H3K9me3 (green), 
H3K27me3 (green) and mH2Al (red), and those that are underrepre- 
sented; H3K4me2 (green) and H3K9Ac (green), (c) Examples of ele- 
vated levels of BRCAl, Y-H2A.X and SNF2H at the Xi (red). SNF2H 
was never observed enriched at the Xi 5/\Z/8 -/-. 



WSTF is one of several proteins that transiently asso- 
ciate with the Xi during S -phase [15]. Other transient 
features of the Xi during this time include the associ- 
ation of the breast cancer gene (BRCAl) and the phos- 
phorylation of histone H2A.X at serine 139 (y-H2A.X) 



[45]. We sought to determine if these features were al- 
tered in the BAZIB knockout cells. Towards this aim, we 
employed IF using antibodies to BRCAl, y-H2A.X, and 
SNF2H— WSTFs partner in the WICH complex [16]. We 
found that BRCAl and y-H2A.X still localize to the Xi 
(1.5% and 2.5% respectively, n = 200) (Figure 2c), similar 
to what is reported in the parental cells [15,45]. In con- 
trast, SNF2H was no longer detected at elevated levels at 
the Xi (n = 1242). However, SNF2H protein still exhibited 
diffuse staining at other locations within the nucleus, con- 
sistent with its role in other complexes independent of 
WSTF and WICH [46-49]. 

Cells lacking WSTF exhibit a unique nuclear phenotype 

DNA staining with 4',6-diamidino-2-phenylindole (DAPI) 
revealed a nuclear phenotype unique to BAZIB knockout 
cells. In parental hTERT-RPEl nuclei, the most intense 
DAPI staining structure corresponds to the Barr body, as 
indicated by the white arrow in Figure 3a. In contrast, the 
nuclei of BAZIB knockout cells show additional DAPI 
dense bodies of comparable to greater intensity as demon- 
strated by the pixel saturation in Figure 3b and 3c. Not all 
nuclei show the appearance of these structures. We ob- 
served that 23.21% of nuclei (n = 2356) possessed numer- 
ous large D API-dense regions (Figure 3d). Intriguingly, 
transient reduction of either component of the WICH 
complex by RNAi was reported to result in a gain of het- 
erochromatin [18], consistent with these observations. 

We sought to further characterize both the compos- 
ition and dynamics of the DAPI- dense regions, as the 
appearance of blocks in only a subset of cells raised the 
possibility that the blocks could be linked to a cell cycle, 
arrest in proliferation, or potentially early stages of cell 
death. 

DAPI-dense regions are enriched for constitutive 
heterochromatin markers 

Areas of intense DAPI staining in nuclei typically corres- 
pond to regions of heterochromatin, such as the consti- 
tutive heterochromatin of mouse major satellite [50], 
senescence associated heterochromatin blocks [51], or 
facultative heterochromatin of the Xi [14]. The blocks 
seen in the BAZIB homozygous mutants were reminis- 
cent of the heterochromatic Barr body in both DAPI in- 
tensity and size. For this reason, we examined the blocks 
for association with heterochromatic features and vari- 
ous complexes known to interact with heterochromatin. 
To accomplish this, we pursued several avenues of in- 
vestigation. First, we explored the relationship of the 
blocks with the three different HPl proteins: HP la, 
HPlp, and HPly [52-54]. Chromatin associated HPl 
levels were shown previously to increase in cells with de- 
pleted levels of WSTF, particularly HPlp [18]. Consist- 
ent with this observation, we found that all three HPl 
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Figure 3 BAZ1B homozygous knockout cells exhibit a unique nuclear phenotype. Representative images of RPEl (a), BAZ1B-H6 (b) and 
BAZlB-04 (c) nuclei counterstained witli DAPI. Wliite arrow-lieads in (b) and (c) point to blocl<-containing nuclei. White arrows point to the DAPI- 
dense Barr body. Images were captured using settings defined for RPEl that result in pixel saturation for block-containing BAZIB -I- nuclei, (d) RPEl 
and B/\l]B -I- nuclei stained with DAPI showing the appearance of large DAPI-dense blocks in a fraction of mutant cells, as indicated to the left. 



proteins localize to the DAPI-dense blocks (Figure 4a). 
To validate that all three isoforms did indeed associate 
with the blocks and that this was not a by-product of 
antibody cross-reactivity, we transiently transfected 
EAZIB knockout cells with HPla, HPlp, or HPly green 
fluorescent protein (GFP) expression constructs. We 
found that all three exogenous HPl types, with HPla to 
a lesser extent, localized to the blocks (Figure 4b). 

Next, we continued to characterize the DAPI-dense 
blocks in the BAZIB knockouts through IF to a variety 
of heterochromatin-associated proteins and covalent his- 
tone modifications. Our results showed that, consistent 
with binding of HPl to H3K9me3 [55-57], the blocks 
were clearly defined by this modification (Figure 4c). 

Given the size of the blocks, we wondered if this rep- 
resented an overall increase in H3K9me3 levels. His- 
tones were extracted from parental and WSTF mutant 
cells and H3K9me3 levels detected by Western blotting. 
Despite the apparent increase by IF, H3K9me3 levels 
were comparable between mutants and parental cells 



(Figure 4d). However, given that only 23% of cells show 
blocks, it is possible that increases are masked by the 
77% of normal cells. 

Recently, the structural maintenance of heterochroma- 
tin protein, SMCHDl, was linked to HPl and H3K9me3 
[58]. In agreement with this, SMCHDl associated with 
the DAPI blocks (Figure 4e). With the exception of the 
Xi, comparable large blocks of HPla, HPlp, HPly, 
H3K9me3 and SMCHDl were not a feature of parental 
cells (Additional file 1). 

Although they did not completely colocalize, areas of 
MBDl and HMG-I/Y enrichment were present within 
portions of the blocks (Additional file 1). As expected, 
the blocks had no association with euchromatin markers 
(Figure 4f), However, 14 other heterochromatin proteins 
did not associate with the blocks, including TIF IB, 
EZH2, Mi2, RbAp46, RbAp48, SAP18, SAP30, SIN3A, 
SIN3B, HDACl, HDAC2, MBD3, ATRX, and SETDBl, 
nor histone HI or the histone variant macroH2Al 
(Additional file 1). 
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Figure 4 BAZ1B -/- DAPI-dense blocks are characterized by constitutive heterochromatin marlcers. (a) Association of endogenous HPla, 
HPip, and HPly (green) witli DAPI-dense blocl<s. Nuclei are counterstained with DAPI (blue) and merged images are shown in the right-side 
panel, (b) Association of exogenous GFP-tagged HPla, HP1|3, and HPly (green) with DAPI-dense blocks, (c) Examples of elevated levels of 
H3K9me3 (green) at DAPI-dense blocks, (d) Western blot showing comparable levels of H3K9me3 in mutant compared to RPEl. Anti-histone H4 
is used as a loading control, (e) Example of the association of SMCHDl (green) with DAPI-dense blocks, (f) Examples showing the lack of associ- 
ation of euchromatin markers H3K4me2 and H3K9Ac with the DAPI-dense blocks. 



DAPI-dense blocks do not correspond to common nuclear 
structures 

We sought to determine if the DAPI-dense regions were 
associated with disruption of B-WICH activity, due to its 
participation in ribosomal DNA transcription at nucleoli 
[20]. The nucleolar organizing region (NOR) of the acro- 
centric chromosomes (13, 14, 15, 21, and 22) localizes to 
and can be used to define the extent of nucleoli [15]. 
We performed DNA FISH using a BAG from the NOR 



as a probe, and IF using an antibody raised against 
MYBBPIA, a member of the B-WICH complex [20], to 
examine if the DAPI blocks coincided with nucleoli. We 
found that neither NOR or MYBBPIA overlapped with 
the DAPI-dense blocks (Figure 5a). 

Given the large number of DAPI-dense blocks, we 
wondered if they corresponded to the inappropriate 
spread of centromeric heterochromatin. The reasoning 
behind this assumption is that cells treated with RNAi 
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Figure 5 BAZ1B -I- DAPI-dense block relationship to other nu- 
clear structures, (a) DAPI-dense block distribution relative to nucle- 
oli markers MYBBP1A (red, IF) and NOR (red, direct-labeled FISH), (b) 
DAPI-dense block distribution relative to centromere markers CENP- 
A (green, IF) and alpha-satellite DNA indicated by a pan-centromeric 
direct-labeled FISH probe (green), (c) DAPI-dense block distribution 
relative to PML bodies (green, IF). 



to BAZIB showed an increase in HPlp, which was pos- 
tulated to represent gain of heterochromatin [18], and 
that HPl is a feature of centromeric heterochromatin 
[59]. We addressed this two ways, first by DNA FISH 
with a pan-centromeric probe, and second by IF to the 
centromeric protein CENP-A [60,61]. We found that the 
blocks did not overlap with either centromeric DNA or 
CENP-A-positive regions (Figure 5b), although signals 
often resided at the edge of the blocks. 

Prior studies have also linked WSTF function to sev- 
eral cellular processes including DNA replication 



(WICH), transcription (B-WICH), and epigenetic silen- 
cing (WICH) [62]. Like WSTF, promyelocytic leukemia 
(PML) bodies have also been linked to the processes of 
DNA replication, transcription, and epigenetic silencing 
[63]. In addition, the number, and sometimes size, of 
PML bodies are reminiscent of the DAPI-dense blocks 
that we see in BAZIB knockout cells. Altogether, this 
led us to question whether the DAPI-dense blocks were 
expanded PML bodies. To test this question, we de- 
tected PML bodies by IF. We found that PML bodies 
had not increased in size and did not colocalize with the 
DAPI-dense blocks (Figure 5c). 

DAPI-dense blocks in BAZIB knockout cells are not linked 
to and do not disrupt the cell cycle 

The appearance of the DAPI-dense blocks in only a sub- 
set of cells could be explained by their appearance being 
associated with a particular period during the cell cycle. 
In order to test this, we examined cells labeled with 
markers of the cell cycle: CDTl for Gl [64,65], PCNA 
distribution for S -phase [66,67], and phosphorylation of 
histone H3 at serine 10 (HSSlOph) for G2 phase [68]. 
DAPI-dense blocks were detected at all stages of the cell 
cycle (Figure 6a). 

Next, we wondered if the appearance of these blocks 
affected the progression of the cell cycle. However, cell 
cycle profiles determined by fluorescence activated cell 
sorting (FACS) were indistinguishable between parental 
RPEl and BAZIB knockout cells (Figure 6b). 

These observations suggest that perhaps cells were 
slowly acquiring blocks over time, or that only a subset 
of cells always have blocks and others do not. 

Cells containing DAPI-dense blocks do not accumulate 
over time in culture and are not stably inherited 

In order to determine if the number of cells with DAPI- 
dense blocks was accumulating over time, we monitored 
cells for the percentage of DAPI-dense positive nuclei 
every 10 passage doublings for a total of 70 doublings. 
We found that at any given passage doubling, the popu- 
lation of cells with DAPI-dense blocks remained consist- 
ently around 20-25% of cells (Table 2). 

Since only a subset of cells possessed the DAPI blocks, 
we questioned whether they were a stable heritable trait 
that could be passed on to daughter cells. If this were 
the case, we would expect that clones derived from a 
non DAPI-dense block containing cell would show no 
blocks, while those derived from a DAPI-dense block 
containing cell would all contain blocks. To investigate 
this question, we isolated single cell clones and exam- 
ined eight subclones of each of three independent 
BAZIB homozygous knockout mutants for the presence 
or absence of the DAPI-dense blocks. Interestingly, all 
subclones contained cells with and without DAPI-dense 
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Figure 6 DAPI-dense blocks in BAZ1B -I- are present in G1, S and G2 phase and do not disrupt the cell cycle, (a) Examples of DAPI-dense 
blocks in BAZIB -I- cells at different stages of the cell cycle as indicated by the appropriate markers; CDTl for Gl (green, top), PCNA for S-phase (red, 
middle) and HSSlOph for G2 (green, bottom). HSSlOph accumulates first at heterochromatic regions before spreading throughout the nucleus [68], 
accounting for the elevated levels observed at the heterochromatin blocks in the example shown. The percentages given in brackets to the left of the DAPI 
image indicate the number of cells positive for the cell cycle stage marker in the population of cells examined (CTDl , n = 83; PCNA, n = 1 05; H3S1 Oph, 
n = 1 15). Note that the numbers are greater than 100%. This is because some H3S10ph positive cells are in late S-phase in RPEl [15], and CDTl, while 
present throughout the cell cycle, is most abundant during Gl [64]. Consequently some cells from other stages of the cell cycle are inadvertently included 
in these numbers, (b) Representative FACS profile for RPEl and 8/\Z/8 -/-. 



Table 2 Appearance of DAPI-dense blocks does not change 
with passage doubling 



Passage 
doubling 


Percentage of nuclei 
with DAPI-dense blocks 


Number of 
nuclei scored 


14 


31.85 


135 


24 


22.74 


277 


34 


22.94 


327 


44 


20.27 


291 


54 


21.30 


315 


64 


24.46 


372 


74 


22.14 


262 


84 


24.13 


377 



Table shows the percentage of BAZIB -/- clone D5 nuclei that showed DAPI- 
dense blocks over the course of 70 passage doublings. A total of 2356 nuclei 
were scored. The mean and standard deviation is 23.73% +/-3.56. 



blocks (Table 3), which suggested that cells could de- 
velop and resolve blocks regardless of the state of the 
parent cell However, it remains a possibility that those 
cells acquiring blocks arrest in the cell cycle and eventu- 
ally die. 

DAPI-dense blocks are dynamic 

The transient transfection of HPlp into BAZIB knockout 
cells clearly highlighted the DAPI-dense blocks (Figure 4b), 



Table 3 DAPI-dense block incidence in single cell clones 



BAZ1B 


Percentage of nuclei with DAPI-dense 


Mean 8f standard 


-/- 


blocks for 8 independent subclones 


deviation 


clone 






D5 


78, 50, 76, 78, 44, 58, 40, 48 


59.0% +/- 16.0 


F3 


8, 40, 22, 56, 24, 42, 20, 24 


29.5% +/- 14.0 


Ml 


20, 28, 20, 18, 22, 20, 32, 38 


24.8% +/- 7.0 



Table shows the percentage of nuclei with DAPI-dense blocks for each of eight 
single-cell subclones derived from three independent BAZIB -/- clones. Fifty 
nuclei were scored for each subclone. 
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providing an opportunity to monitor live cells expressing 
the GFP construct in real time. Cells expressing HP1|3- 
GFP were followed over 24 hour and 72 hour periods. 
Since the number of block-containing cells did not in- 
crease in the population over time, the blocks could either 
be spontaneously forming and resolving, or those cells 
with blocks could arrest, or undergo apoptosis. We ob- 
served that the block-containing nuclei did not undergo 
apoptosis, but were readily able to divide, ruling out arrest 
in the cell cycle (See Additional file 2). We found that a 
block-containing cell was able to give rise to a daughter 
cell that also acquired blocks, and, therefore, the condition 
was mitotically heritable (Additional file 3). Furthermore, 
the blocks appear to resolve with the exception of the Xi- 
associated HP1|3-GFP block, highlighting the dynamic na- 
ture of the blocks. Interestingly, we noted that, despite 
substantial movement by the cell body, the blocks were 
static within nuclei (Additional file 4), similar to observa- 
tions made in real-time with replication factories [67]. Fi- 
nally, blocks were observed for as long as 24 hours in cells 
before dividing, and appeared to resolve right before enter- 
ing mitosis, confirming that they are not transient features 
of a particular cell cycle stage (Additional file 3). Given 
these observations and the fact that single cell clones give 
rise to a similar fraction of cells with DAPI-dense blocks 
as the parental cell population (and not 100% block- 
containing or 100% block-free), we conclude that the ap- 
pearance and resolving of the blocks is a dynamic process. 
What triggers the formation and subsequent resolution of 
the blocks in some cells and not others is presently 
unknown. 

Introduction of exogenous full length WSTF ORF Is 
sufficient to remove the appearance of DAPI-dense blocks 

We wondered if the DAPI-dense blocks in BAZIB 
knockout cells could be removed by reintroducing wild 
type WSTF protein. To test this, an expression construct 
containing a Myc-tagged full-length WSTF ORF was in- 
troduced into WSTF knockout cells by Nucleofection. 
Cells were examined by IF using antibodies specific to 
the Myc-tag and WSTF. Cells that were positive for Myc 
were also positive for WSTF (n = 94) (Figure 7a). Fur- 
thermore, transfected cells were restored for transient 
association of WSTF to the Xi as defined by H3K27me3 
[15] (Figure 7b). We examined Myc-positive nuclei for 
DAPI-dense blocks. All Myc-positive nuclei showed 
DAPI patterns comparable to parental RPEl, with no 
evidence of intense DAPI-dense blocks (n=119), sug- 
gesting that the blocks can be resolved by reintroduction 
of wild-type protein. 

Loss of WSTF results in gene expression changes 

Given that the DAPI-dense blocks are enriched for 
markers of constitutive heterochromatin, and that WSTF 
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Figure 7 Distribution of exogenous Myc-tagged WSTF in BAZIB 
-/- nuclei, (a) Distribution of exogenous WSTF as detected by anti- 
WSTF IF in BAZIB -I- nuclei (green) (top row) and the MYC-tag (red) 
(bottom row). Examples of elevated levels at the Xi are indicated with 
the white arrowheads, (b) Elevated levels of exogenous WSTF-Myc 
(red) at the Xi as defined by H3K27me3 (green). The location of the Xi 
is indicated by the white arrowhead. Nuclei are counterstained with 
DAPI (blue). 

k ) 

is involved in transcription [62], we assessed the WSTF 
mutant cells for changes in gene expression, with the 
thought that perhaps most gene expression changes 
would be reductions due to the gain in heterochromatin. 
Total RNA was extracted from three separate cultures of 
parental RPEl cells as well as from three separate cul- 
tures for each of three independent BAZIB knockout 
clones, and one heterozygous mutant clone. Residual 
DNA in the RNA preparation was digested and its re- 
moval confirmed by PGR before preparing cDNA, label- 
ing and hybridizing to a NimbleGen human gene 
expression 12xl35K array containing >45,000 transcript 
targets. Genes that had significant changes of two fold 
or more compared to RPEl (p < 0.05) were identified for 
each mutant. The heat map shown in Figure 8a, clearly 
shows consistent and substantial changes in gene ex- 
pression between the three independent BAZIB mutants 
as compared to the parental RPEl cells. The genes that 
were not shared between all three independent mutants 
were excluded. The remaining genes of interest can be 
found in Additional file 5 (n = 515). In the knockout 
cells, genes for which expression levels changed were 
not all reductions; some went up (173) while others 
went down (342). To validate the microarray data, we 
performed quantitative RT-PGR (qRT-PGR) on several 
genes that showed increases or decreases in the micro- 
array data relative to the parental RPEl cells. All 
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Figure 8 Gene expression changes in BAZ1B -I- cells, (a) 

Dendrogram (Heat Map) of hierarchical clustering of gene expression 
between RPEl and three independent BAZIB -I- clones. Heat map 
represents 586 data points, corresponding to probes showing greater 
than 2-fold change in gene expression at 95% confidence. Expression 
levels range from +1.4 to +14.2, corresponding to signal intensity on the 
microarray and represented by a blue (low) to red (high) gradation indi- 
cated by the scale bar beneath the heat map. Each column represents 
the averaged signal from three replicate data sets for each clone. Each 
row represents a single specific probe. Changes in gene expression can 
be seen by the difference in row colors between the mutants and paren- 
tal cells, (b) Representative examples of genes activated in 5/\Z/8 -/- 
cells by standard end-point RT-PCR analysis shown as an inverted image 
of an ethidium bromide stained agarose gel. The gene name is given to 
the left, and the cell line indicated at the top. Lanes show with (+RT) and 
without (-RT) reverse transcriptase and a water control (H20). 



consistently showed the same trend (data not shown). 
Interestingly, of the 515 gene expression changes ob- 
served in the BAZIB knockouts, 72 of these were shared 
in the heterozygous mutant cells (Additional file 6). This 



indicates a degree of sensitivity to WSTF haploinsuffi- 
ciency, further supporting an important contributing 
role oi BAZIB hemizygosity in WBS. 

Despite our supposition that most genes would be si- 
lenced due to the gain in heterochromatin, as mentioned 
above, a substantial number of genes showed increased 
expression in the mutants relative to the parental cells. 
Among these, a subset appeared to show activation in 
both BAZIB -I- and BAZIB +/- cells (Figure 8b). 
These genes could not be assessed by qRT-PCR, as fold 
change cannot be determined when expression of a con- 
trol or unknown sample is zero, but the endpoint RT- 
PCR analysis is consistent with the microarray data. 

Gene expression changes were not restricted to certain 
chromosomes or chromosomal regions (Additional file 5). 
With the exception of chromosome 21, all chromo- 
somes have genes that showed increases or decreases in 
expression, although all changes were less than 2.0% of 
the total gene content for the chromosome (coding and 
non-coding). Given that chromosome 21 has the lowest 
gene content, it does not seem that this chromosome is 
an exception for any biological reason. Despite the 
blocks frequently residing close to nucleoli, the NOR- 
containing chromosomes (13, 14, 15, 21 and 22) did not 
have the highest number of silenced genes. Instead, 
the chromosomes that contained the largest number of 
silenced genes were chromosomes 20 (1.55%) and 12 
(1.50%). No obvious trend was observed for silenced 
genes concerning genomic location, such as residing 
close to large satellite blocks or in pericentromeric or 
subtelomeric regions. Interestingly, the protocadherin 
cluster on chromosome 5q31.3 that is associated with 
various neuropsychiatric disorders [69], and implicated in 
tumor suppressor activity [70,71], consistently showed an 
increase in expression of most of the alpha members, but 
decreases in the adjacent beta and gamma clusters. Intri- 
guingly, in the heterozygous BAZIB mutant, the alpha 
cluster was up-regulated, similar to the complete knock- 
out, but changes to the beta and gamma clusters were not 
observed. 

Discussion 

Gene silencing at the Xi is remarkably stable, involving 
multiple redundant gene silencing mechanisms that en- 
sure silencing is maintained throughout all subsequent 
somatic cell divisions. Stability of the imprinted Xi in ex- 
traembryonic material is compromised upon loss of the 
Polycomb group protein Embryonic Ectoderm Develop- 
ment (EED) [72,73], but there are few other examples of 
Xi instability due to a single gene loss. In fact, global re- 
activation of the Xi has only been achieved through fu- 
sion of female mammalian cells with mouse embryonal 
carcinoma cells [74-77] or via reprogramming mouse fe- 
male cells to pluripotency through forced expression of 
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transcription factors [78] . Therefore, if WSTF is involved 
in maintaining chromatin states at the Xi, it is not sur- 
prising that the loss of WSTF did not result in an obvi- 
ous wave of X-linked gene reactivation and loss of 
bonafide Xi markers, such as XIST or H3K27me3. Genes 
that are reactivated from the Xi are unlikely to be 
expressed at levels dramatically different from their Xa 
counterpart. Consequently, the anticipated change in 
gene expression would likely be less than two fold. How- 
ever, reassessing the data for 1.5 fold change did not 
show a bias toward increase in X-linked gene expression 
levels compared to the autosomes (data not shown). Al- 
ternatively, however, these observations can be inter- 
preted to suggest that perhaps the transient association 
of WSTF with the chromosome, as the Xi DNA is repli- 
cated [15], is unrelated to maintenance of XCI. 

WSTF, as part of the WICH complex, has been shown 
to transiently associate with the Xi prior to y-H2A.X and 
BRCAl [15]. The fact that BRCAl and y-H2A.X continue 
to transiently associate with the Xi when WSTF is absent 
suggests that the two cell-cycle linked events are not 
coupled. Indeed, we have previously speculated that 
WICH may transiently associate with the newly replicated 
DNA in order to reestablish the basal state of H2A.X- 
Y142ph. This theory supports the idea that WSTF recruit- 
ment to the Xi is independent of y-H2A.X and BRCAl as- 
sociation since the ground state of H2A.X-Y142ph is set 
independently of DNA damage [21]. 

Intriguingly, almost one quarter of BAZIB knockout 
cells displayed numerous, large, DAPI-dense blocks. Con- 
sistent with the DAPI staining, the blocks are character- 
ized by constitutive heterochromatin-associated proteins, 
including HPl, HMG-I/Y, and MBDl, and by the covalent 
histone modification H3K9me3. In addition, the HPl- 
interacting protein SMCHDl [58], was also seen localizing 
to the blocks. While the blocks did not overlap with nucle- 
oli, they were often seen abutting these nuclear structures, 
as indicated by NOR DNA FISH and MYBBPIA IF; 
markers of nucleoli [15]. Prior studies have shown that, in 
addition to rDNA-containing regions of the acrocentric 
chromosomes, heterochromatin domains of most human 
chromosomes (including regions found at telomeres, cen- 
tromeres, and chromosome arms), associated with nucle- 
oli [79]. Perhaps, perinucleolar regions are sites of 
heterochromatin maintenance. Indeed, in mouse, the 
binding partner of Wstf, Snf2h, transiently accumulates in 
a perinucleolar region as the mouse Xi is replicated [80]. 
Perhaps, in BAZIB knockout cells, checks on heterochro- 
matin formation are compromised resulting in the in- 
appropriate spread into genomic regions not normally 
associated with heterochromatin. The obvious gain in the 
large blocks of heterochromatin support a model in which 
WSTF assists in maintaining chromatin states, and in its 
absence, too much heterochromatin forms [17]. 



The size and number of DAPI-dense blocks in BAZIB 
knockout cells indicates that the gain of heterochroma- 
tin in affected cells is substantial. Previously, it was 
shown that artificially reducing WSTF protein levels by 
RNA interference, was associated with lower levels of 
transcription, as measured by monitoring the levels of 
bromouridine incorporation into RNA [18]. By micro- 
array analysis, a large number of genes consistently 
showed a reduction in expression (342) in three inde- 
pendent mutants. This change could be attributed to the 
gain in heterochromatin, although secondary effects, 
such as silencing of a key transcription factor, cannot be 
ruled out. Of interest, the reduced expression was not 
global, impacting only a subset of all genes. This sug- 
gests that if WSTF is acting to reestablish appropriate 
chromatin states post DNA replication, as postulated 
[17], it is either only acting at specific regions of the gen- 
ome, as most genes are unaffected by the proteins loss, 
or that certain regions of the genome are more vulner- 
able to inappropriate heterochromatin spread/assembly 
than others. 

We also observed a large number of genes that 
showed an increase in expression (173) in all three 
BAZIB knockout clones. This is consistent with WSTF s 
previously established role as a regulator of transcrip- 
tion, capable of activating or silencing genes [62]. Not- 
ably, loss of a single BAZIB allele was sufficient to 
impact a subset of those genes affected in the knockouts, 
indicating that some genes are more sensitive to WSTF 
levels than others. This has important implications for 
BAZIB hemizygosity in WBS, and our heterozygous tar- 
geted clones provide an important model in which to 
pursue the contribution of WSTF to the phenotype in 
this complex contiguous gene deletion disorder. 

Some genes showed reactivation from a silent state in 
the mutants. In contrast, relatively few genes showed the 
opposite trend of complete silencing in the mutants. It is 
possible that within the pool of cells assessed, complete 
silencing might only be observed in cells with DAPI- 
dense blocks that represent fewer than 25% of the total 
cell population. In these cells, genes embedded in the 
blocks or within their vicinity might be fully silenced; 
yet not in the 75% of cells that appear normal. Continu- 
ing with this notion, it is possible that the actual number 
of genes impacted by WSTF loss might be greater than 
our microarray data would otherwise suggest. If some 
genes show increases or decreases in gene expression 
that are close to, but greater than, our 2-fold cutoff only 
in DAPI-dense block nuclei, these changes might be re- 
duced to below the cutoff by a lack of change in the ma- 
jority non-DAPI-dense block cells. 

Despite the extensive size of the DAPI blocks, we did 
not see an increase in nuclear levels of H3K9me3. 
Therefore, an alternative explanation for the reactivation 
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of genes could be an indirect effect, as H3K9me3 is 
redistributed to the blocks, essentially reducing local re- 
pressive effects elsewhere. 

The dynamic nature of the blocks is intriguing. Cells 
that acquire blocks do not undergo apoptosis and the 
cell cycle is apparently unaffected. Single cell cloning in- 
dicates that block formation is reversible, and that the 
overall number of cells with blocks remains relatively 
constant within an actively dividing population and be- 
tween independent mutants. Live cell imaging confirms 
that block appearance is not linked to the cell cycle, and 
indicates that blocks are essentially immobile, and do 
not prevent cells from proceeding through mitosis. Re- 
introduction of wild-type WSTF resolves the blocks, and 
therefore provides a powerful manipulable system in 
which to investigate regions of the WSTF protein that 
are required for this activity. 

Conclusions 

In this study, we have explored the impact of WSTF 
loss. Our findings emphasize that many aspects of the 
function of WSTF remain enigmatic. Identif)^ing the 
common characteristics of the DNA sequences in the 
DAPI-dense blocks may provide clues as to how WSTF 
contributes to gene regulation. It is important for future 
studies to examine the mechanism of WSTF action so 
that we may fully understand the role that WSTF plays 
in nuclear processes and its contribution to WBS. In 
addition, our heterozygous mutant cells provide a useful 
model to investigate the effects of BAZIB hemizgosity 
and WSTF haploinsufficiency in the absence of the loss 
of the other 27 genes seen in WBS. 

Methods 

Cell culture 

The cell line used to generate the knockout clones was 
hTERT-RPEl, a 46, XX telomerase-immortalized human ret- 
inal pigment epithelial cell line (Clontech, Laboratories, Inc., 
Mountain View, CA, USA: No. C4000-1). The cells were 
maintained according to the supplier s recommendations. 

Generating BAZIB knockout cell lines 

The targeting template for homology mediated repair 
was generated as follows. The left homology arm (HA) 
was amplified by PGR with the primers LHA-F2 and 
LHA-R2 (Additional file 7) using RPEl genomic DNA as 
a template. The same was done for the right HA with 
the primers RHA-Fl and RHA-Rl (Additional file 7). 
The PGR products for the left and right HAs were TA 
cloned into pDRIVE (Qiagen, Valencia, GA, USA). Plas- 
mid DNA was harvested using the Nucleospin Plasmid 
kit (Machery-Nagel, Bethlehem, PA, USA), before se- 
quence verifying. The insert of the left HA was excised 
from the pDRIVE vector using Notl and Xbal, and the 



right HA was excised using Xhol and Mlul. The pSEPT 
cassette was digested with Xbal and Xhol and the pBG- 
KO-tk-A vector (Ghadwick, manuscript in preparation) 
was digested with Notl and Mlul, All digests were sepa- 
rated on agarose gels (Agarose Unlimited, FL, USA), and 
inserts excised and purified using a Nucleospin Gel and 
PGR Glean-up kit (Machery-Nagel, Bethlehem, PA, 
USA). The left HA, pSEPT cassette, and right HA, were 
ligated overnight at 16°G into pBG-KO-tk-A vector using 
T4 DNA ligase (NEB, Ipswich, MA, USA). The ligation 
was transformed into NEB- 5 alpha chemically competent 
E. coli cells according to the supplier s recommendations 
(NEB, Ipswich, MA, USA). Plasmid DNA was prepared 
using the Plasmid Plus Midi kit (Qiagen, Valencia, GA, 
USA). The plasmid DNA was linearized by digestion 
with Ascl restriction enzyme, and gel purified. To gener- 
ate BAZIB knockout cell lines, we introduced ZFN ex- 
pression constructs (Sigma- Aldrich, St. Louis, MO, USA: 
Gat No. GSTZFN-IKT) and the linearized repair tem- 
plate by Nucleofection (Lonza Group Ltd, Basel, 
Switzerland). Twenty-four hours post-Nucleofection, 
cells were seeded into four 96-well plates in media sup- 
plemented with neomycin (150 (ig/mL). After approxi- 
mately two weeks, clones were expanded and genomic 
DNA isolated using a Wizard SV 96 Genomic DNA Puri- 
fication System kit (Pr omega, Madison, WI, USA). Gor- 
rectly targeted clones were screened on the left and right 
sides by PGR between a primer to the genomic locus, but 
distal/proximal to the HA, and a primer located in the 
pSEPT cassette. The screening pairs were Left-Screen-F2 
with pSEPT-Rl and pSEPT F2 with Right-Screen-Rl 
(Additional file 7). Glones that generated the correct sized 
PGR product using these two screening pairs were tested 
for NHEJ by PGR using the primers BAZIB Sigma-F and 
BAZIB Sigma-R (Additional file 7). If a product was gen- 
erated, it was sequenced to determine if NHEJ had re- 
sulted in insertions or deletions. All restriction enzymes 
were obtained from NEB (Ipswich, MA, USA). 

WSTF rescue 

A full-length wild-type human myc-tagged WSTF expres- 
sion construct was obtained from Origene (Rockville, MD, 
USA: Gat.No. RG216159). The expression construct was 
introduced into BAZIB knockout cells by Nucleofection 
using a 4D-Nucleofector (Lonza Group Ltd, Basel, 
Switzerland) and solution SF on program EN- 138 accord- 
ing to the manufacturers recommendations. 

Antibodies 
Immunofluorescence 

Rabbit anti-WSTF antibodies were obtained from Gell 
Signaling Technology (Danvers, MA, USA: Gat No. 
2152) and Bethyl Laboratories (Montgomery, TX, USA: 
Gat No. A300-446A). Rabbit anti-SMGHDl (Gat No. 
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A302-872A) was also purchased from Bethyl Laborator- 
ies (Montgomery, TX, USA). Rabbit anti-CDTl (Cat No. 
10/2012), rabbit anti-SETDBl (Cat No. 06/2012), rabbit 
anti-CENP-A (Cat No. 2186S), mouse anti-TIFlB (Cat 
No. 5868S), and mouse monoclonal antibodies to myc 
(Cat No. 2276S) were purchased from Cell Signaling 
Technologies (Danvers, MA, USA). Mouse monoclonal 
antibodies BRCAl (Cat No. sc-6954) were obtained from 
Santa Cruz Biotech (Dallas, Texas, USA) as were goat 
anti-SAP18 (Cat No. sc-8473), mouse anti-Histone HI 
(sc-8030), goat anti-MBDl (Cat No. sc-9395), goat anti- 
MBD2 (Cat No. sc-12444), goat anti-MBD3 (Cat No. 
sc-9402), goat anti-RbAp46 (Cat No. sc-8272), goat anti- 
RbAp48 (Cat No. sc-8270), goat anti-ATRX (Cat No. 
sc-10080), rabbit anti-SIN3A (Cat No. sc-994), rabbit 
anti-SIN3B (Cat No. sc-768), rabbit anti-Mi2 (Cat No. 
sc-1564), and rabbit anti-HMG-I/Y (Cat No. sc-1564). 
Antibodies to Mouse monoclonal antibodies to 
phospho-H2A.X (Serine-139) (Cat No. 05-636) and 
SNF2H (Cat No. 05-698) were obtained from Millipore 
(Billerica, MA, USA), as were rabbit polyclonal anti- 
bodies to EZH2 (Cat No. 07-400), histone H3 trimethy- 
lated at lysine-27 (Cat No. 07-449), histone H3 
phosphorylated at serine 10 (Cat No. 06-570), histone 
H3 dimethylated at lysine 4 (Cat No. 07-0030), histone 
H3 acetylated at lysine 9 (Cat No. 06-942), and histone 
H3 trimethylated at lysine 9 (Cat No. 07-523). Mouse 
anti-HPla (Cat No. MAB3584), anti-HPl|3 (Cat No. 
MAB3448), and anti-HPly (Cat No. MAB3450) were 
obtained from Chemicon International (EMD Millipore, 
Billerica, MA, USA). Mouse anti-PML (Cat No. 
ab96051) antibody was purchased from Abeam (Cam- 
bridge, MA, USA). Mouse anti-MYBBPlA (Cat No. 
SAB1400390) was obtained from SIGMA-Aldrich (St. 
Louis, MO, USA). Rabbit anti-macroH2A.l antibody 
was described previously [81]. Alexa-Fluor® conjugated 
secondary antibodies were obtained from Life Technolo- 
gies Corporation (Invitrogen, Grand Island, NY, USA). 

DNA and RNA FISH probes 

Direct-labeled FISH probes were prepared using a Nick 
Translation kit with SpectrumGreen dUTP according to 
the manufacturer s recommendations (Abbott Molecular, 
Abbott Park, IL, USA). PI -artificial chromosome (PAC) 
RP5-1174A5, which contains rDNA and detects the Nu- 
cleolar Organizing Regions (NOR), was obtained from 
Dr. P. Finelli (Finelli et al. 2012). The pan-centromeric 
probe (Cat No. 1695-F-02) was purchased from Cambio 
(Cambridge, UK). The XIST RNA-FISH probe was as 
described previously [15]. 

Western blotting 

The protein samples for Western blotting were prepared 
by lysing approximately one million cells in RIPA buffer 



according to standard procedures [82] . Samples were sep- 
arated on 4-20% acrylamide gels at 200 V. Histone sam- 
ples for Western blotting were prepared according to 
Abcams histone extraction protocol (http://www.abcam. 
com/index.html?pageconfig=resource&rid=l 1410). Sam- 
ples were separated on AnykD acrylamide gels (Biorad, 
Hercules, CA, USA) at 200 V. All SDS-PAGE gels were 
transferred to pre-equilibrated 0.2 [i PVDF (Biorad, 
Hercules, CA, USA) using a BioRad Turbo Transfer ap- 
paratus. The PVDF membranes were blocked with 5% 
non-fat dry milk in 0.1% PBS-Tween 20 for 30 minutes 
and washed twice in 0.1% PBS-Tween 20/1% igepal for 
15 minutes each. Primary antibodies were diluted in 
blocking solution and either incubated with gentle rocking 
overnight at 4C (WSTF Cell Signaling Technologies and 
Bethyl Laboratories antibodies) or for 1 hour at room 
temperature (all other antibodies). The membranes were 
washed twice as described above. The membranes were 
incubated with horseradish peroxidase-conjugated (HRP) 
secondary antibodies diluted in blocking solution for 
1 hour at room temperature with gentle rocking. The 
membranes were washed as previously described. Detec- 
tion of HRP was achieved by incubating membranes with 
equal volumes of ECL Solution A and Solution B (GE 
Healthcare, Pittsburgh, PA, USA) for 3 minutes at room 
temperature. The two WSTF antibodies described above 
were used. Rabbit anti-B-Actin (Cat No. 11/2012) was 
purchased from Cell Signaling Technology (Danvers, MA, 
USA). The rabbit anti-H3K9me3 antibody described above 
for IF analysis was used for Western analysis. Rabbit anti- 
Histone 4 (Cat No. 07-108) was obtained from Millipore 
(Billerica, MA, USA). HRP anti-rabbit secondary antibody 
was obtained from Cell Signaling Technology (Danvers, 
MA, USA: Cat No. 7074S). 

DNA and RNA FISH and indirect immunofluorescence 

To prepare slides for IF and FISH studies, approximately 
2x10^ cells were grown directly on standard microscope 
slides. IF was performed as described [15]. DNA was coun- 
terstained using ProLong® Gold antifade reagent sup- 
plemented with 4', 6-diamidino-2-phenylindole (DAPI) (Invi- 
trogen, Grand Island, NY, USA). For IF combined with DNA 
FISH, IF was first performed. The FISH probes were dena- 
tured in the presence of human Cot-1 DNA (Invitrogen, 
Grand Island, NY, USA) at 78°C for 10 minutes followed by 
30 minutes at 37°C to block repetitive sequences before ap- 
plying to the slides. Slides were then treated as previously de- 
scribed [15]. Antifade containing DAPI was added to the 
slides, followed by a cover slip. IF combined with RNA FISH 
was performed essentially as described [83]. 

Image acquisition 

Images were collected using a DeltaVision pDV (GE 
Healthcare, Pittsburgh, PA, USA). Delta Vision images 
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were deconvolved with softWoRx 3.7.0 (Applied Preci- 
sion, Issaquah, WA, USA) and compiled with Adobe 
Photoshop CS2 (Adobe Systems, San Jose, CA, USA). 
Live cell imaging videos were recorded using an Andor 
Revolution Spinning Disk Laser Confocal microscope 
and an Andor 897 EMCCD camera (Andor Technology 
pic. Belfast, UK). Live cell imaging data was processed 
using Andor iQ2 software. 

HP1-GFP Expression constructs 

GFP-tagged HPl-a, |3 and y expression constructs were 
generated by excision with EcoRl of the full ORF inserts 
from previously described pcDNA3.1-CT-Myc con- 
structs [84], into the EcoRl site of pcDNA3.1-CT-GFP 
(Invitrogen, Grand Island, NY, USA). 

Fluorescence activated cell sorting (FACS) analysis 

For FACS analysis, an asynchronously growing popula- 
tion of 1 X 10^ cells were stained with propidium iodide 
[85] and flow sorted using a FACSCanto (BD Bioscience, 
San Jose, CA, USA). Cell cycle profiles were generated 
using FACSCanto software. 

Microarray gene expression analysis 

Microarray expression analysis was performed on paren- 
tal cells as well as three independent knockout mutants 
and one heterozygously targeted clone. For each sample, 
total RNA was extracted from three independent cell 
pellets isolated at different time points. Total RNA was 
extracted from cells using the RNeasy Mini Kit from 
Qiagen (Valencia, CA, USA). Fluted RNA was treated 
with DNasel (NEB, Ipswich, MA, USA) before cleaning 
and purifying the RNA using an RNeasy Mini Kit col- 
umn (Qiagen, Valencia, CA, USA) according to the man- 
ufacturers recommendations. RNA quality was assessed 
by agarose gel electrophoresis and absence of contamin- 
ating DNA was confirmed by PCR with primers to the 
large tandem repeat DXZ4 [86]. Single stranded cDNA 
was prepared using the High Capacity cDNA Reverse 
Transcription Kit (Invitrogen, Grand Island, NY, USA: 
Cat. No. 4368814). The cDNA was then provided to 
the FSU-Nimblegen Microarray Facility for labeling, 
hybridization, and washing. Probes were hybridized onto 
a Human Expression 12xl35K array (Roche-Nimblegen, 
Indianapolis, IN, USA: 05543789001). Hybridization data 
was analyzed using DNA ArrayStar version 4 (DNAS- 
TAR, Madison, WI, USA). The three replicates of RPEl 
were compared to one another. Transcript targets were 
included if their confidence level was > 95% based on a 
p-value of p < 0.05. This was done for each of the mu- 
tant clones. Replicate sets of each of the mutants were 
compared to the RPEl replicate set. In order to be in- 
cluded, greater than two fold change in expression was 
required, with a confidence level of >95% based on a 



p-value of p < 0.05. Additionally, transcript targets that 
did not meet these criteria in all three of the independ- 
ent mutants were excluded. Microarray expression data 
has been submitted to the Gene Expression Omnibus, 
and has been assigned accession number GSE47390. 

Quantitative RT-PCR Analysis 

All qRT-PCR was performed on a CFX96 (Biorad, 
Hercules, CA, USA) using EvaGreen 2X qPCR Master- 
mix (ABM, Richmond, BC, Canada: Cat No. Mastermix- 
S). The two-step amplification cycle began with a 10 mi- 
nute incubation at 95C, followed by 40 cycles of 15 sec- 
onds at 95C and 30 seconds at 60C. Fluorescence was 
measured at the end of each of the 40 cycles. After- 
wards, a melt curve was generated by measuring the 
fluorescence after heating for 5 seconds in 0.5C incre- 
ments, beginning at 65C and ending at 95C. qRT-PCR 
analysis of BAZIB used Quantitect primer assay 
QT00072044 along with GAPDH assay QTOl 192646. 
Both were obtained from Qiagen (Valencia, CA, USA). 

Standard RT-PCR Analysis 

Standard RT-PCR was performed using OneTaq (NEB, 
Ipswich, MA, USA). The cDNA samples were denatured 
at 94°C for 30 seconds, amplified through 40 cycles of 
30 seconds at 94°C, 30 seconds at 58°C, and 30 seconds 
at 68°C, and finally held at 15°C. Details on the primers 
used can be found in Additional file 7. 

Additional files 



Additional file 1: Representative IF images of chromatin proteins 
relative to the DAPI-dense blocks and the normal distribution of 
H3K9me3, HP1 and SMCHDI in parental RPEl cells. Top panels show 
representative IF images of BAZIB -I- cells with DAPI-dense blocks, 
showing the distribution of various chromatin proteins (red or green). Bot- 
tom panels show the normal distribution of H3K9me3, HPl and SMCHDI in 
parental RPEl nuclei. The nucleus is counterstained with DAPI (blue). 

Additional file 2: Movie clip of a HPip-GFP Nucleofected BM'IB -I- 
heterochromatin block containing cell entering mitosis uninhibited. 

Movie shows a cluster of three BAZIB -I- clone D5 nuclei that are 
positive for HP1(3-GFP blocks. Toward the end of the movie, the middle 
cell enters and exits mitosis. The movie represents a time lapse of 
21 hours and 50 minutes with images captured every 5 minutes under 
40X magnification. Note: The duration of the cell cycle is potentially slo- 
wed due to frequent exposure to the illumination source. 

Additional file 3: Movie clip of a HPip-GFP Nucleofected BAZ^B -I- 
heterochromatin block containing cell entering mitosis and giving 
rise to a block-containing daughter cell. Movie shows a BAZIB -/- 
clone D5 nuclei that is positive for HP1(3-GFP blocks (top left), that passes 
through mitosis and gives rise to a block-containing nucleus. The blocks 
appear to resolve toward the end of the movie leaving the Xi-associated 
HP1(3-GFP as the only remaining obvious block. For comparison, an 
HPl (3-GFP positive non-block containing nucleus shows only the Xi- 
associated HPl (3-GFP (lower central portion of movie). The movie represents 
a time lapse of 62 hours and 20 minutes with images captured every 10 
minutes under 60X magnification. Note: The duration of the cell cycle is 
potentially slowed due to frequent exposure to the illumination source. 
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Additional file 4: Movie clip of a HPip-GFP Nucleofected BAZ1B -/- 
heterochromatin block containing cell showing the static nature of 
the block positions in the nucleus. Movie shows a BAZIB -I- clone 
D5 nuclei that is positive for HPlp-GFP blocks, that develops intense 
blocks that show no obvious movement despite the mobility of the cell. 
The movie represents a time lapse of 62 hours and 20 minutes with images 
captured every 10 minutes under 60X magnification. Note: The duration of 
the cell cycle is potentially slowed due to frequent exposure to the 
illumination source. 

Additional file 5: Excel file listing genes that consistently change 
expression by >2-fold in BAZIB -I- cells. File contains tabs for all data 
combined, as well as gene expression for each chromosome individually. 

Additional file 6: Excel file listing genes that consistently change 
expression by >2-fold in BAZIB +/- cells that are shared with 
BAZIB -I- cells. 

Additional file 7: Table listing all oligonucleotides used in this study. 
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tyrosine 142; H3K4me2: Histone H3 dimethylated at lysine 4; H3K9Ac: Histone 
H3 acetylated at lysine 9; H3K9me3: Histone H3 trimethylated at lysine 9; 
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